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We investigate cooling and trapping of single atoms inside an optical cavity using a quasi-resonant 
field and a far-off resonant mode of the Laguerre-Gauss type. The far-off resonant doughnut mode 
provides an efficient trapping in the case when it shifts the atomic internal ground and excited 
state in the same way, which is particularly useful for quantum information applications of cavity 
quantum electrodynamics (QED) systems. Long trapping times can be achieved, as shown by full 
3-D simulations of the quasi-classical motion inside the resonator. 
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I. INTRODUCTION 



Cavity quantum electrodynamics (QED) is a powerful tool for the deterministic control of atom-photon interactions 
at the quantum level. In fact, the strong confinement allows to achieve the strong coupling regime where single quanta 
can profoundly affect the atom-cavity dynamics [jjj . Many experiments have now reached this regime and interesting 
phenomena have been demonstrated, as a quantum phase gate JE|, the Fock state generation in a cavity field ||, 
and quantum nondemolition detection of a single cavity photon 0]. Moreover, cavity QED offers unique advantages 
for quantum communication and quantum information processing applications. In fact, atoms may act as quantum 
memories, while photons are flexible transporters of quantum information, and quantum networks of multiple atom- 
cavity systems linked by optical interconnects have been already discussed in the literature || . The primary technical 
challenge on the road toward these applications is to trap individual neutral atoms within a high-finesse cavity for 
a reasonably long time. Recent experiments [^|J^] have already succeeded in trapping single atoms inside an optical 
cavity driven at the few photon level, just using the strong coupling with the cavity QED mode for both cooling and 
trapping. However, the scheme of Refs. HQ is not entirely suitable for quantum communication purposes because 
it has limited operation flexibility and provides short trapping times (order of hundreds of microseconds). In fact, it 
employs a single cavity mode, while it is preferable to have an additional trapping mechanism which does not interfere 
with the cavity QED interactions able to provide the atom-photon entanglement needed for the manipulation of 
quantum information. With these respect, another experiment has already demonstrated significant trapping times 
(~ 28 ms) of single Cs atoms within a cavity, employing an additional far-off resonance trapping (FORT) mode [||. 
O^' Using the fact that, in the strong coupling regime, the trajectory of an individual atom can be monitored in real time 
by the quasi- resonant cavity QED field ||, the FORT beam can be turned on as soon as the atom enters the cavity 
■ ■ *h . in order to increase its trapping time. 

Several mechanisms for cooling inside an optical resonator have been already discussed in the literature Il0|jll] , p^ , [l3| , 
involving either cavity mode driving, or direct atom driving via a classical laser field from the side, or even active 
feedback on atomic motion as in [l4| . However, only the recent paper by van Enk et al. has discussed in detail the 
effects of an additional FORT beam on the cooling and trapping dynamics and its interplay with the quasi-resonant 
cavity QED field. Here we shall consider a situation analogous to that of Ref. |p| , even though we shall extend our 
study to new FORT mode configurations. We have chosen a parameter region corresponding to the weak driving 
limit, with an empty-cavity mean photon number N e ~ 0.01. We have performed full 3-D numerical simulations of 
the quasi-classical atomic motion, including the effects of spontaneous emission and dipole-force fluctuations. 

In two-level systems, red-detuned FORT beams shift the atomic excited state |e) up and the ground state \g) 
down by the same quantity, and this is the most common situation, studied in great detail in Jl5[ . However, the 
most interesting situation for quantum information processing applications is when both levels are shifted down by 
the FORT beam: in this case, excited and ground state atoms are trapped in the same position, and this greatly 
simplifies the quantum manipulation of the internal state. In fact, the most flexible situation for quantum information 
processing is having a trapping mechanism independent of the atomic internal state. This configuration can be realized 
by using a FORT that is red-detuned in such a way that the excited state is relatively closer to resonance with a 
higher-lying excited state than with the ground state. We shall discuss in detail both situations (equal or opposite 
optical Stark shifts), and we shall find the interesting result that in the case of equal Stark shifts for ground and 
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excited state, the use of doughnut modes, that is, higher-order Gauss-Laguerre modes, as FORT mode, is able to 
increase significantly the trapping time within the cavity. 

The paper is organized as follows. In Section II we describe the physics of an atom trapped in an optical potential 
and strongly interacting with a cavity QED field. We also discuss the changes introduced by the FORT doughnut 
mode. In Section III, applying the general approach developed in Ref. |ll| to a cavity mode configuration, we discuss 
the conditions under which the center-of-mass motion of the atom can be adiabatically separated from the internal and 
cavity mode dynamics, and treated in a quasi-classical way. The corresponding 3-D Fokker-Planck equation for the 
phase-space atomic motion will be derived. In Section IV, the results of the numerical simulations of the corresponding 
stochastic differential equations will be presented in detail, and in Section V these results will be discussed for both 
equal and opposite energy shifts of the |e) and \g) states. Section VI is for concluding remarks. 



II. THE PHYSICAL PROBLEM 



We consider a two level atom coupled to a quantized cavity mode and to an additional classical red-detuned FORT 
beam, coinciding with another longitudinal mode of the cavity, with a wavelength As longer than that of the quasi- 
resonant cavity mode, X g . The common situation is to consider lowest order Gaussian modes for both fields [fl5|| , 
having their maximum intensity along the cavity axis. Here we shall consider a different situation, where the FORT 
mode is a higher-order Gauss-Laguerre mode, the so-called doughnut mode, having its maximum intensity at a nonzero 
radial distance from the cavity axis. This means that the atoms are trapped out of the cavity axis (see Fig. |l| for 
a schematic description of the system). At first sight, this choice may look not optimal, because in this case the 
coupling with the fundamental Gaussian cavity QED field responsible for cooling is smaller. Nonetheless, we shall see 
that this choice is convenient in the case when the classical FORT mode shifts down both excited and ground state, 
which is the most interesting case for quantum information processing applications. 




FIG. 1. Schematic description of the system. Atoms falling from a magneto-optical trap are trapped within a cavity with a 
fundamental Gaussian quantized resonant mode (yellow) and an intense classical doughnut FORT mode (light blue). 



We recall that the Laguerre-Gauss modes LG pm are the solutions of the paraxial Helmoltz equation in cylindrical 
coordinates (p, t9, x) 1 17|. In this paper we consider the doughnut modes with radial index m ^ and azimuthal index 
p = 0, whose intensity is given by 
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where P is the power, k$ = 2ir /Xg the wavenumber, and W{x) the beam radius. The doughnut mode radius is given 
by the position of the radial maximum, given by p max = W(x)y / m/2. We can simplify the description assuming a 
nearly planar cavity so that W(x) ~ Wg. The classical, red-detuned FORT mode induces AC Stark shifts AE g on 
the ground state \g) and AE e on the excited state |e) We can distinguish two different situations: (a) \g) is 

shifted down and |e) is shifted up, which happens in the common situation where only the e g transition is close to 
the red detuned FORT mode. In this case the two shifts are opposite, (b) Both \g) and |e) are shifted down, which 
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happens when |e) is closer to resonance with a higher-lying excited level than with the level \g) pjj , p(| , pl| . We 
shall study both situations and in case (b) we shall assume that the various detunings can be chosen so that the two 
levels are shifted down by the same quantity. Furthermore, to simplify the comparison between the two cases, we 
shall assume that the two situations occur with the same wavelength As. 

Considering a frame rotating at the probe driving frequency tv p , the Hamiltonian of the system can be written as 

IPI 2 

H = l - + hujap^a + hwgptfa + V(R) + h (Ea) + E*a) , (2) 

where uj ap = uj eg — uj p is the detuning of the the atomic resonance from the probe frequency, uj gp = uj g — oj p is the 
detuning of the cavity QED mode with annihilation operator a, cr = \g)(e\, E is the cavity driving rate, and i?, P 
are the position and momentum vector operators of the atom, having a mass M. The interaction potential V(R) 
describes the interaction between the internal atomic levels, the cavity modes and the atomic center-of-mass motion, 
and it is given by the coupling with the quantized cavity mode and with the FORT doughnut mode. Making the 
usual dipole and rotating wave approximations, this interaction term can be written 

V(R) = Tig(R) (oV + a^a) + V s (R), (3) 

where 

g(p,x) = 50 sinOgx) exp j-^jj (4) 

is the space-dependent Rabi frequency due to the coupling with the quantized mode, and the second term V (R) 
describes the effect of the Stark shifts induced by the FORT mode, assuming the following form in the two cases, (a) 
and (b): 

case (a) V s (R) = fiS{R){^a - aa^) (5) 
case (b) V s {R) = -hS(R). (6) 



The frequency ac-Stark shift is generally given by S(p, x) = aIo m (p, x)/4h 22 23y, where a is the atomic polarizability 
and Iom(p, x) is the FORT intensity of Eq. ([!]), so that one can write 

S(p, x) = S 0P 2m sin 2 (k s x) exp j - ^ J , (7) 

with So = 2 m+1 P a / TThWg m+2 . We have generally considered different waists W g and Ws for the two modes, but 
they practically coincide in typical situations. 

The FORT mode provides the main trapping mechanism. This means that the atom will be trapped around the 
anti-nodes of the red-detuned FORT field, because So > go and the cavity QED field is weakly driven. However, 
due to the different wavelengths, one does not have a periodic situation, and the atom feels a different cavity QED 
coupling in different wells. In the experiment of Ref. || the cavity length L is such that 2L — 104A 5 = IO2A5. To 
simplify our simulation we have however chosen L — 16A g = 15As, as it has been done also in |pl| . This is equivalent 
to choose a fictitious larger value for A^ which however does not modify the essential physics of the problem, because 
the FORT mode is in any case far-off resonance. With this choice we consider only 30 potential wells in the cavity, 
only 8 of which are quantitatively different (see Ref. [^5[ ). 

Dissipation, diffusion and all non-conservative effects appear due to spontaneous emission and cavity losses. The 
quantum evolution of the atom-cavity system is therefore described by a master-equation for the atom-cavity density 
operator p 0] 

= — — [H, p] + k (2apa) — a) ap — pa) a) — 7 (a^ap + pa^a) 

+ ?1 [d 2 k D(k ■ x)e-^ n ap^e^ R (8) 

47T J 

where n is the cavity damping rate and 7 is spontaneous emission decay rate. The last term describes the effects of 
atomic recoil, with k giving the direction of emitted photon and D(k ■ x) = (1 + (k ■ x) 2 )/2 describing the angular 
pattern of dipole radiation p3 . 
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III. QUASI-CLASSICAL DESCRIPTION OF THE ATOMIC MOTION 



The internal and cavity dynamics is governed by the detunings, oj ap and oj gp , the spontaneous emission rate 7, 
the driving rate E, the cavity damping rate k, the FORT shift So, and the coupling constant go- In the parameter 
region relevant for current experiments, this dynamics is much faster than the atomic motional dynamics, especially 
for heavy atoms as Cs || or Rb [Q . Therefore, internal and cavity dynamics can be adiabatically eliminated in order 
to obtain a reduced effective evolution equation for the motional degrees of freedom only. This adiabatic elimination 
can be made in a systematic way by adapting the general approach developed for free space cooling situations in 
Ref. [jl6| to the present cavity scheme (see also fOjfl). One starts from the evolution equation of the Wigner operator 
obtained by performing the Wigner transform only on the motional Hilbert space, and writing this equation as a 
Taylor expansion in terms of two small parameters of the problem. Denoting with k ~ k g ~ ks — k a = oj eg /c the 
typical wavenumber of the problem, we have that one small parameter is 

showing that the motional state is characterized by a momentum spread Ap much larger than the momentum kicks 
felt by the atom during any photon emission or absorption. The second small parameter is 

which instead shows that the Doppler shift associated with the momentum spread is small with respect to the atomic 
and cavity linewidths. The two conditions set a lower and an upper bound for the momentum spread of the atom, 
which, thanks to the Heisenberg inequality, put also a lower and an upper bound for its position spread. These bounds 
allow to describe the atomic motion in classical terms, because the atom is sufficiently localized in phase space to 
make it possible to describe its motion in terms of trajectories, while at the same time the phase space spread always 
remains larger than h (see [p|Jll|]). The crucial point is that the two conditions (|^) and [l0|) must be consistent and 
this happens when 

-C Try, Kk. (11) 

This means supplementing the well known necessary condition for laser cooling in free space, h 2 k 2 /2M <C fvy (the 
atom has to be still in resonance after spontaneously emitting) |HJ, with an analogous condition for the exchange of 
cavity photons. 

For small parameters t\ and £2, the atomic motion is much slower than the internal one, and the adiabatic elimination 
of the latter is obtained if the atomic kinetic energy term is neglected in the equation for the Wigner operator, which 
then effectively factorizes into a product of a stationary state r\ (f) for the cavity and internal degrees of freedom 
evaluated at the fixed atomic position r = (p,6,x), and a real- valued motional Wigner function f(r,p,t). The 
stationary state ij (r) satisfies the steady-state Bloch equation 



LBioch(r)v('r) = ~i 



u! nn a' er + cu nn a' a 



:lP Ja + (Ed + E*a) + -p-,r){r) 



+ k (2ar)(r)a ii — arj(r) — rj(r)a a) 

+ 7 (2cr?7(f)cr t - £rVj?(r) - r](r)a^a) = 0, (12) 

where LBioch(r) comes from a zeroth-order expansion of the master equation (|g) in the small parameters ei and £2 of 
Eqs. (|) and ©. 

The resulting equation for the Wigner function f(r,p,t) is of Fokker-Planck type and is given by fl3| , |l6|| 

i 1 ij 3 

+ E D ^)g^f(r,P> t)+J2 T ^)^-(P 3 mP, t)). (13) 
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All the coefficients of this Fokker-Planck equation depend upon average values and correlation functions of the internal 
and cavity degrees of freedom, evaluated on the stationary state at fixed atomic position rj (r) of Eq. ( |l2| ) . An example 
is provided by the average value (cV) (r) appearing in Eq. ( |l3| ) in the diffusion term due to spontaneous emission, 
which also depends upon the diagonal matrix En given by E xx = 2/5, E yy — E zz — 3/10 [^3lp3. The other atom- 
cavity quantities determining the Fokker-Planck equation coefficients are $ = a^a + a^a and ^ = a^a — au^ , whose 
expectation value determines the mean dipole force acting on the atom </>(r). In fact, 

^_/-^(r)($)(r)-^5(rlW(r) case (a) 

n> I -h§g{r)^){r) + hf F S{r) case (b) ' 1 > 

which is a sum of the quantized mode contribution and the FORT mode contribution. This is no more true for the 
friction matrix I\j(r) and the diffusion matrix Dij(r), which are bilinear functionals of the dipole force operator and 
assume different forms in the two cases (a) and (b). In fact, in case (6) of equal shifts, the FORT mode simply adds 
a conservative potential, independent of the internal atomic state, giving no contribution to friction and diffusion 
(provided that the FORT beam intensity fluctuations are negligible, see p2|,p3||). Therefore we can write 



r??(r-Hrf/(r)+rf/(r) case (a) 
if? (r) case (b) 



D 



, W _J^+<W + flg S (f1 case(a) . . 



D??(r) case(b) 



where the quantized mode contribution (gg), the FORT mode contribution (SS), and the cross contribution (gS) 
have been singled out. These terms can be then written as 



and 



D « (f> = s^w^wrw (20) 



(21) 



(o = ^(oj^o^o + h2 ^^ cS9{ff) 

^(0 = ^(0^(^(0, (22) 



where 



X 39(r)=j/ drr([*(r),$(0)]) (23) 
Jo 

/>OC 

X gS (?)=i drr([$(T),*(0)]) (24) 
Jo 

/>oo 

X S9 (r)=i dTr([*(r),$(0)]) (25) 
Jo 

/>oo 

x ss {r f )=i dTT([*(r),*(0)]), (26) 



and 



£"(0 = / rfr 



o 



o 



1({<I>(t),$(0)})-(<J>} 2 



(27) 
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£, Ss (r) = 



dr 



dr 



dr 



({3(t),¥(0)}>-<$> <*) 



2<{*(r),$(0)})-<*) ($) 



2<{*(r),*(0)})- W 



(28) 
(29) 
(30) 



The cross diffusion term proportional to Tfoj (r) is usually much smaller than the diffusion terms due to the dipole force 



fluctuations Dij(r) and due to spontaneous emission oc Eij |13| , |16| , and we shall neglect it. Using this approximation 
and the above expression for the diffusion matrix, the six-dimensional phase space diffusion matrix of the Fokker- 
Planck equation (H|) becomes semipositive-definite. This means that it can be associated to a classical phase space 
stochastic process, describing the stochastic trajectories of the atomic center-of-mass within the cavity. Consistently 
with the adiabatic and quasi-classical description discussed above, the motional Wigner function f(r,p,t) becomes 
therefore a nonnegative, classical phase space probability distribution [ fl6| |. Our numerical analysis is based just on 
the simulation of these stochastic 3-D trajectories, which are obtained as solutions of the Ito stochastic equations 
associated to the Fokker-Planck equation (fi~3j). Moreover, as it can be easily seen from Eq. (|), (@), and Eqs. @- 



(|22|), one has T xx oc \ g • 



T zz oc Wa , and the same is true for the 



diffusion matrix Dij. Since it is always X g , As <C W g , Ws, it is evident that the only relevant term in the friction and 
diffusion matrices is the xx component along the cavity axis, where both the quantized field and the FORT mode 
show the largest spatial gradients JD|. In our numerical simulations, we have therefore considered both the friction 
force and the dipole force contribution to diffusion, only along the cavity axis x, while we have kept the spontaneous 
emission diffusion terms in all three directions. 

Taking into account Eq. (13), and the above approximations, we have therefore numerically solved the following 
Ito-equations p5|,p6[|: 



dr 



P_ 

M" 



dt 



dp = <f)(r)dt — 



+ fik a \jl (<rV) (r) 




VD^ATjdWi 





(31) 



where dW\ and dW = (dW x ,dW y ,dW z ) are four independent, zero mean, Wiener increments with the property 
dWidWj — 2Sijdt, = 1, x y z). The quantities $(r), T xx (r), D xx (r), (crV) (f) in Eq. ( ^H ) have been determined 



by numerically solving Eq. (|12j) in the atom-cavity Hilbert space truncated at n — 4 photons (see J27J for details). 
We have also checked that these numerical solutions reproduce the results of the analytical approach of [0 in the 
weak-driving limit. 



IV. NUMERICAL RESULTS 



In the next subsections we present and discuss the results of the numerical simulations in both cases, (a) and 
(b). We have chosen equivalent conditions for the two cases, so that the corresponding numerical results are directly 
comparable. 



A. Study of case (a) 



We consider parameter values referred to the experiment of Ref. gj. In fact, we consider the Cs transition at 
A a = 2n/k a — 852.4 nm between the ground state \g) = \6Si/ 2 ,F = 4,mp — 4) and the excited state |e) = 
\§P?,/2;F — h.rap — 5), and a cavity mode resonant with it, i.e., uj ap — uj gp — — A p = 2ir x 10 MHz (A p is the 
probe detuning from resonance). The spontaneous emission rate is 7 = 2ir x 2.6 MHz, while the other quantized 
mode parameters are g — 2tt x 30y/e MHz, n = 2n x 4 MHz, W g = 20 /mi, and E = 6.77 MHz, so that the empty 
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cavity mean photon number is N e — E 2 / (k 2 + A 2 ) = 0.01. Since it is go > 7, k, we are therefore in the strong 
coupling regime of cavity QED. Then we use a LG01 doughnut mode as red-detuned FORT field, with parameter 
values Ws = 20 /im, Sq — ne/2 MHz//im 2 . These choices give p max ~ 14.1 pm for the doughnut radius, and the 
maximum Stark shift given by the FORT mode, achieved at p max , is S'max = 27r x 50 MHz. 

As discussed above and in p5[ , we have chosen a fictitious large value of As, so that L = 15As = 16A S (X g = X a ) in 
order to simplify the simulation, without however changing the physics because the exact value of As is unimportant 
as long as it is far-off resonance. The dominant potential is the one due to the FORT and therefore the atomic 

equilibrium positions will be situated at the FORT antinodes x e n = (n — l/2)As/2, n — 1, , along the axial 

direction and at the nonzero radial distance from the cavity axis p max ~ 14.1 /im. Because of the sharp radial 
potential well of the doughnut mode, the atom will experience a sharp radial confinement around p max besides the 
axial confinement around an antinode. 

The superposition of the two modes with different wavelengths determines a spatially aperiodic situation within 
the cavity. The above choice gives 8 not equivalent potential wells in which the atom is subject to different couplings 
with the quantized cavity mode. This aperiodic situation is shown in Figs. § and ^, where the axial friction coefficient 
r xx (r) of Eqs. ( ^B| ) and (|l7|)-(|l^), and the dipole contribution to the axial diffusion coefficient D xx (r) of Eqs. (|l^ ) 
and (E^)-(]22), at fixed radial distance p = p max (bottom of the radial well), and as a function of the axial coordinate 
x, are plotted. 




-1.5 - 



-2 ' • • • • • • • 1 

0.5 1 1.5 2 2.5 3 3.5 4 

x/A.s 

FIG. 2. Axial friction coefficient T xx (r), at fixed radial distance p = /5 max (bottom of the radial well), versus the rescaled 
axial coordinate x/Xs, in the case (a) of opposite Stark shifts induced by the FORT mode. Parameter values are in the text. 




FIG. 3. Dipole contribution to the axial velocity diffusion coefficient D xx (r)/M 2 , at fixed radial distance p — /9 ma x (bottom 
of the radial well), versus the rescaled axial coordinate x/Xs, in the case (a) of opposite Stark shifts induced by the FORT 
mode. Parameter values are in the text. 
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1. Simulation results 



To characterize the trapping and cooling dynamics, we have carried out a quantitative study inside a central well, the 
one centered in x% — 2.25As and ranging from 2As to 2.5A5. In order to simulate the typical experimental condition, 
we have chosen proper initial conditions, considering that, in the experimental procedure, the FORT field is switched 
on only when the laser probe transmission exceeds a fixed threshold indicating the presence of the atom inside the 
cavity. The initial position has been taken axially A5/8 away from equilibrium point x\ (then xq — 2.125As), radially 
along the doughnut maximum po = p max , and uniformly distributed over the polar angle 0. For what concerns the 
initial velocity, it is reasonable to choose a vertical velocity with components v x q = v y o = 0, v z q — 10 cm/s. 

In order to examine qualitatively the typical atomic motion, we report some snapshots from two simulated trajec- 
tories with slightly different initial conditions: one describes a tangential incidence of the atom along the doughnut 
perimeter (xq = 2.125As, yo = Pmax, z o = 0, v x q = v v q = 0, v z q = 10 cm/s) (see Figs. (§|6|); the other instead describes 
an orthogonal incidence with respect to the doughnut mode (xq — 2.125As, yo — 0, zq — p max , v x q = v v q = 0, v z q = 10 
cm/s) (see Fig. [?]). 
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FIG. 4. Time evolution of the axial position in the case of an initial velocity tangential with respect to the doughnut FORT 
mode. Parameter values are in the text. 
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FIG. 5. Time evolution of the radial coordinate in the case of an initial velocity tangential with respect to the doughnut 
FORT mode. Parameter values are in the text. 
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FIG. 6. Radial trajectory in the case of an initial velocity tangential with respect to the doughnut FORT mode. Parameter 
values are in the text. 
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FIG. 7. Radial trajectory in the case of an initial velocity orthogonal with respect to the doughnut FORT mode. Parameter 
values are in the text. 



The trajectory is well confined, axially around a FORT antinode, and radially around the maximum /O max : the 3D 
atomic motion occurs substantially on a plane orthogonal to the cavity axis, with the trajectory drawing just the 
shape of the doughnut mode. From Figs. ^-^| one can recognize that the atomic motion is characterized by three 
different time scales: 

1. The fastest timescale is given by the axial oscillations, which have, for our parameter values, a time period 
- 2^s (see Fig. |). 

2. A slower timescale is associated with the radial oscillations, characterized by a time period ~ 100 fis (see Fig. ||); 
these oscillations become wider in the case of orthogonal incidence with respect to the doughnut mode, since, 
in this case, the atom probes more the doughnut radial elasticity (compare, in fact, Fig. || where the oscillation 
amplitude is ~ 1 /im with Fig. [?], where the amplitude is ~ 4 /im). 

3. The third and slowest timescale is given by the atomic rotations around the cavity axis. For tangential incidence, 
the initial angular momentum is large and the rotation period is ~ 1 ms while, for orthogonal incidence, the 
atom acquires a nonzero angular momentum only because of radial diffusion, and the rotation period is larger, 
~ 10 ms. 



The final escape of the atom from the cavity is practically always along the axial direction, and this is due to the 
heating provided by the strong axial diffusion, which prevails with respect to the radial diffusion, determined only by 
the spontaneous emission. 
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To determine the mean trapping times, we have defined as trapping time T the time spent by the atom inside 
a single potential well, As/2 wide along the axial direction, and with a radius equal to 2Ws- Sampling about 400 
simulated trajectories, we have found the results displayed in Figs, g and 0. 
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10 1 10 2 10 3 10* 10 5 

trapping time (|is) 

FIG. 8. Root mean square velocity along the cavity axis v™ 3 as a function of the trapping time T, for each simulated 
trajectory, in the case (a) of opposite Stark shifts. 

In Fig. U we have displayed the root mean square velocity along the cavity axis v™ s as a function of the trapping 
time T for each simulated trajectory. We can see a clear separation of simulated points. In fact, about 60% of atoms 
are not trapped at all: they correspond to the upper set of points in Fig. 0, having a velocity v™ s > 20 cm/s, and 
for which the trapping time is below 2 ms. These are not cooled via the cavity QED interaction, and the velocity 
y rms ^ 28 cm/s is mainly a result due to the initial conditions. These uncooled atoms are those more influenced by 
the spatial region where the axial friction is negative (see Fig. 0), where atoms can be accelerated. The remaining 
40% of atoms (the points below the threshold of 20 cm/s in Fig. ]8]) are trapped: their velocity v™ s of about 14 cm/s 
is a result of the cooling provided by the exchange of cavity photons. These atoms reach thermal equilibrium and 
have trapping times greater than 2 ms. 

Considering only the subset of trapped atoms, they have a probability P(t) to be trapped for a time greater than 
t. The trapping time statistics is shown in Fig. [)[ where it is also compared with a decaying fitting curve P th (full 
line in Fig. ||). The best fitting mean trapping time is r = J* °° t (—dP th /dt) dt = (17 ± 1) ms, which is comparable 
to the experimental value obtained using a fundamental Gaussian FORT mode in Ref. []3| and with the numerical 
simulations performed in |lq| , again for a fundamental Gaussian FORT mode. This is not surprising because, except 
for the fact that the atom is now trapped at a nonzero distance from the cavity axis, the physics of cooling is similar 
to that occurring in a lowest-order Gaussian FORT mode, and the analysis of Ref. |l5| can be essentially repeated. 



0.4 




2 4 6 8 10 

time (us) x10 * 

FIG. 9. Trapping time statistics for the subset of trapped atoms in the case (a) of opposite Stark shifts. P(t) is the 
probability for an atom to be trapped for a time greater than t. The full line is a fitting decay curve yielding a mean trapping 
time r = (17 ± 1) ms. 
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B. Study of case (b) 



We consider the same parameter values of case (a), except that we slightly adapt the probe detuning and choose a 
value A p = —2ir x 35 MHz. As a consequence, we then set E = 22.13 MHz, in order to keep the same empty-cavity 
mean photon number N e = 0.01 of case (a). 

The situation is in many respects very similar to the preceding one: the atomic equilibrium positions are again 
situated at the FORT antinodes x e n — (n — l/2)As/2, n — 1, ... along the axial direction and at the nonzero radial 
distance from the cavity axis p max ~ 14.1 /im. There is still an aperiodic situation with 8 not equivalent wells within 
the cavity. However, for equal atomic shifts, the FORT does not affect both friction and diffusion (it does not affect 
the force fluctuations, see Eqs. ( |l5| ) and (|l6|)) and therefore, for what concern friction and diffusion, an axially periodic 
situation is restored in this case, with a spatial period set by the cavity QED mode. The periodic friction and diffusion 
spatial variations are shown in Figs. |l^ and O, where the axial friction coefficient T xx (r) of Eqs. ( |l5| ) and (p^7|)-(^9|), 
and the dipole contribution to the axial diffusion coefficient D xx (r) of Eqs. ( |l6| ) and (pp[)-(p2[), at fixed radial distance 
P = Pmax, and as a function of the axial coordinate x, are plotted. 



1000 - 



2 2.5 

xA, s 



FIG. 10. Axial friction coefficient F xx (f), at fixed radial distance p = p max (bottom of the radial well), versus the rescaled 
axial coordinate xj As, in the case (b) of equal Stark shifts induced by the FORT mode. 
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FIG. 11. Dipole contribution to the axial velocity diffusion coefficient D xx (f) /M , at fixed radial distance p = p max (bottom 
of the radial well), versus the rescaled axial coordinate x/Xs, in the case (b) of equal Stark shifts induced by the FORT mode. 



If we compare Figs. ^ and || with Figs. [T(] and |TT] we see that in case (b) the maxima of the axial friction and 
diffusion coefficients are lower, but one has the advantage that now the friction coefficient is always positive, while in 
case (a) it may assume very large negative values, yielding heating rather than cooling. This always positive friction 
implies that all atoms falling in the cavity are now cooled and trapped, while in case (a), a fraction of the atoms can 
be heated and are not trapped. 
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1. Simulation results 



In order to calculate the mean trapping time in case (b), we have carried out a computer simulation considering 
the same initial condition of case (a) (see section [ V A 1 ) . Sampling about 400 simulated trajectories, we have found 
the results shown in Figs. [L2 and O. 



0.3 - 



0.15 - 



0.1 - 
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trapping time (us) 

FIG. 12. Root mean square velocity along the cavity axis v™ a as a function of trapping time T, for each simulated trajectory, 
in the case (b) of equal Stark shifts. 

In Fig. |l2] we have displayed the root mean square velocity along the cavity axis v^ ms as a function of trapping time 
T for each simulated trajectory. At variance with case (a), now all atoms are trapped, with trapping times greater 
than 1 ms. This is due the fact that axial friction is always positive and atoms are cooled everywhere. In this case 
the probability P(t) to be trapped for a time greater than t (shown in Fig. |l3| ) is computed considering all simulated 
points. The data are again fitted by a decaying fitting function P th {t) (full line in Fig. |l^), yielding a mean lifetime 
t = | °° t (-dP th /dt) dt = (44 ± 3) ms. 




2 

time (us) x10 5 

FIG. 13. Probability P(t) for an atom to be trapped for a time greater than t in the case (b) of equal Stark shifts. The full 
fine is a fitting decay curve, yielding a mean trapping time r = (44 ± 3) ms. 



This is the most relevant result of our investigation, showing that using a doughnut mode as FORT mode allows to 
achieve significant trapping times in the case (b), when the FORT induces equal Stark shifts on the atomic levels. In 
fact, we get a mean trapping time larger than that obtained, in the same situation, with a red-detuned fundamental 
Gaussian mode (see the numerical analysis of Ref. fig ], Section IVF, where r = (28 ± 2) ms). This shows that, in the 
case of equal Stark shifts, the radial confinement at a nonzero distance from the cavity axis provided by the doughnut 
FORT mode is preferable with respect to the radial trapping along the cavity axis provided by the TEMrj,o FORT 
mode. The improvement provided by doughnut mode is useful for quantum information processing applications in 
cavity QED systems, because when the FORT mode induces equal shifts on the ground and excited state, the two 
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states are trapped at the same positions within the cavity, and the internal state can be manipulated independently of 
its center-of-mass state. In quantum information applications, it is important to evaluate the variation of the cavity- 
QED coupling Ag. In our parameters' regime it is Ag 16%, which is worser than that for a TEMo,o FORT mode 
under the same conditions (Ag k: 6%). However the problem associated with such a variation of the coupling can be 
circumvented if one adopts an adiabatic transfer scheme such as the one suggested in for quantum information 
processing with hot trapped atoms. In such a case, the operation time T op for the adiabatic transformations must 
satisfy the conditions 2/is ~ < T op -C v~ ad -, ~ 100/is, where v axial, radial are the axial or radial oscillation 

frequencies. 

An interesting alternative to the configuration studied here is to replace the TEMo,o cavity QED mode with the 
higher order LGo.i cavity QED mode, having the same coupling g at p max - In this way we achieve a better overlap 
between the FORT mode and cavity QED one, which leads to a smaller variation of coupling A g , and therefore also 
to a better cooling configuration. In fact, repeating the simulations with the LGo,i cavity QED mode we have seen a 
slight improvement of the trapping time. 

It is clear that the additional practical difficulty of using another LG mode makes the experimental realization even 
harder; however this should be very challenging since, extending the configuration from one LG cavity-QED mode 
to many degenerate LG cavity-QED modes, one has new possibilities, as following atomic motion in detail j30), or 
increasing the cooling effect as shown in [|3l| . 



V. DISCUSSION OF THE RESULTS 



Let us now discuss in detail which are the main features of using red-detuned doughnut modes as FORT fields. 
The atom-cavity dressed picture provides an intuitive way to understand the advantages brought by the use of the 
doughnut mode. In the present case of very weak driving (N e = 0.01 in our case), the cooling mechanism is well 
described in terms of the eigenstates of the atom-cavity system (the dressed states) containing at most one excitation 
(see also pHJl5| , p8t ) . The state with no excitation is the ground state |0), with energy Eq = —hS(r), while the first 
two dressed states with one excitation |±) have energies 



_ J huj a ± fiyjg 2 {r) + S 2 (r) case(a) , . 

± ~ \ huj a - hS(r) ± hg(r) case(b) ' 1 > 

so that the transition frequencies (relative to tu a ) from the ground state |0) to excited states |±) have the expressions 



(E± -Ep) _jS(r)±yJg*(r)+S*(r) case(a) 

±g(r) case(b) 



The spatial variation of these transition frequencies along the axial direction (and at the radial position p = p 
shown in Fig. [H] for the case (a) of opposite Stark shifts, and in Fig. [ll] for the case (b) of equal Stark shifts 



max 
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FIG. 15. Transition frequencies A± (relative to u> a ) from the ground state |0) to the first excited dressed states |±) along the 
axial direction at fixed radial distance p — p ma x (bottom of the radial well), in the case (b) of equal Stark shifts. The straight 
line gives the probe detuning A p . 



The patterns of spatial variation of the transition frequency allow to understand the difference between cases (a) 
and (b). The cooling mechanism provided by the cavity QED interaction in the weak driving limit, can be understood 
in analogy with Doppler cooling. In fact, by tuning below resonance, stimulated absorption of a probe photon followed 
by spontaneous emission or cavity decay leads to a loss of energy. The maximum cooling rate is achieved when the 
excitation rate times the detuning is maximum. If the detuning turns from red to blue, cooling is replaced by heating, 
while if the detuning becomes more red, the atom is still cooled, but at a lower rate. In case (a) (Fig. 14), it is better 
to tune the probe to the lower dressed state because it has smaller spatial variations (see the straight line in Fig. |l4| ) 
and it is easier to reach a compromise between having cavity regions with an optimal cooling rate (red detuning) and 
not too large regions with heating (blue detuning). However, as the atom moves radially from the doughnut intensity 
maximum at p = p max , the situation rapidly worsens, either if the atom moves towards the center (the FORT beam 
decreases while the quantized field increase and the blue detuned region increases) or if tends to leave the cavity (the 
detuning becomes more red and the cooling rate decreases) . 

In case (b), one tunes again to the lowest dressed state in order to have red detuning, and therefore cooling, 
throughout the cavity, when the atom is at the radial equilibrium position p max (see the straight line in Fig. fL5| ). 
From Eq. ( |35| ) we see that the probe detuning decreases (becomes less red) when the atom moves radially towards the 
center, while it becomes more red if the atom moves away radially. However, the probe detuning can be chosen so to 
remain always red: in this way the atom is cooled everywhere, and is never accelerated. For this reason the doughnut 
FORT mode provides longer trapping times in case (b) of equal Stark shifts rather than in case (a) of opposite 
shifts. What is more important is that, in case (b), the doughnut FORT mode provides longer trapping times than a 
fundamental Gaussian FORT mode (see Section IVF). In fact, in the latter case, the probe detuning would be 
tuned so to have optimal cooling on the cavity axis, where the atoms will be now trapped. However, because of radial 
diffusion (due to spontaneous emission) which also leads to an increasing angular momentum and the consequent 
rising of a centrifugal potential, the atom tends to move radially away from the cavity axis, so that the optimal 
cooling condition is rapidly lost. In fact, for increasing p, the axial pattern of Fig. |l5| rapidly vanishes, the driving 
probe becomes too far-off resonance and the cooling efficiency is lost. 

The advantage of the doughnut FORT with respect to the TEMo^o FORT is that it imposes a much sharper radial 
well around the equilibrium position corresponding to the optimum cooling condition. The atom is less free to move 
radially, and moreover is subject to a smaller centrifugal force because it is trapped by the doughnut at a larger radial 
distance from the cavity axis. In other words, the radial potential well provided by the doughnut FORT mode is more 
suitable to counteract the radial departure of the atom and so to preserve the optimal cooling condition. 

The importance of a strong radial confinement to optimize cooling and trapping in the case of equal Stark shifts can 
be illustrated also checking that when the width of the radial well is decreased, the trapping time increases. To this 
purpose, we have carried out a simulation where a higher order Laguerre-Gauss FORT mode, LGo,i2, is used instead 
of the LGo,i FORT mode, keeping the other conditions unchanged. To be more specific, we have considered the same 
experimental parameters of the case with the LGrj.i FORT mode in case (b), in such a way that the LG0.12 mode 
gives the same Stark shift S m!iX = 2tt x 50 MHz at the same radial distance p ma x ~ 14.1 urn. This means adapting 
both the incident FORT power and the doughnut mode waist, so that the only difference with the case studied in 
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the preceding Section is the smaller width of the radial well (see Fig. [T|, where the radial profile of the Stark shift 
induced by the LG0.12 FORT mode is compared with that of the LGo,i mode). The chosen parameters are the same 
as above except that now the spatial dependence of the Stark shift (equal for both ground and excited state) is 



S(p,x) = Sqp sin (ksx) exp 



with W s = 20/V12 ~ 5.77 fim and S ~ 1.2 x 1CT 2U MHz//im 



(34) 




20 25 30 35 40 45 

p im) 

FIG. 16. Stark shifts 5*0,1 and S'0,12 (narrower curve) as a function of the radial coordinate p, at the axial antinode. Parameter 
values are in the text. 



Considering again the initial conditions specified in Section IV A 1 and sampling about 400 simulated trajectories, 
we have found the results shown in Figs. [l7 and [ll| 
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FIG. 17. Root mean square velocity along the cavity axis v™ s as a function of trapping time T, for each simulated trajectory 
in the case (b) and with the higher order FORT mode LGo,i2- 

In Fig. ^ we have displayed the root mean square velocity along the cavity axis v r x ms as a function of trapping 
time T for each simulated trajectory. All atoms are again trapped, with trapping times greater than 1 ms. The 
probability P(t) to be trapped for a time greater than t (shown in Fig. |l^) is computed considering all simulated 
points. The data are again fitted by a decaying fitting function P th {t) (full line in Fig. |r|), yielding a mean lifetime 
t = f™t (-dP th /dt) dt = (64 ± 5) ms. 
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FIG. 18. Probability P(t) for an atom to be trapped for a time greater than t in the case (b) and with the higher order 
FORT mode LGo,i2- The full line is a fitting decay curve yielding a mean trapping time r = (64 ± 5) ms. 

This significant improvement in the mean trapping time, provided by the LGo,i2 FORT mode with respect to the 
LGo^ one, is due to the strongest radial confinement achieved (see Fig. |l6|). The atom stays closer to the radial 
potential minimum where cooling conditions are optimized and the probability to leave the cavity decreases. This is 
a further argument showing that the main advantage of using a doughnut FORT mode instead of a TEM ,o FORT 
mode for atom trapping in the case of equal Stark shifts is just the stronger radial confinement which makes easier to 
optimize the cooling conditions. The more the atom is radially confined, the longer is trapped. 

Here we have decreased the width of the radial well by choosing a higher order doughnut mode and by simultaneously 
adjusting its waist, so to have an unchanged radial equilibrium position p m ax- This solution may be practically difficult 
to implement because one would need different cavity mirrors for this high-order doughnut FORT mode. Anyway, 
even though the implementation of the trapping scheme with a higher-order doughnut mode is difficult, the latter 
numerical results clearly show the importance of the radial confinement, and the direction one has to follow in order 
to increase the trapping time of neutral atoms in cavities in the case of equal Stark shifts induced by the FORT mode. 

A possible solution to overcome the practical difficulties of a LGo,i2 experimental set-up can be the following one, 
provided that one has enough power and stability for the laser driving the FORT mode. The same radial potential 
around p max as that created by the LGo,i2 FORT mode can be achieved using a very intense LGo,i FORT mode. One 
can take So = 47re MHz//jm 2 , leading to S max — 2tt x 400 MHz, and leave all the other parameters unchanged (see 
IV. B for the parameters of the LGo,i FORT mode in case (b)). In this new configuration, both the radial and the axial 
wells are 8 times deeper, and therefore, together with the desired strong radial confinement (substantially the same 
as that of LGo,i2 FORT mode), we have a strong axial confinement. Both these facts lead to a considerable increase 
of trapping time, up to 2 seconds (a result of the same order as the one recently reported in J2f|), as demonstrated 
by the numerical simulation of some trajectories which we show in the following table: 



trapping time r (ms) 


eq. temp. v" ns (cm/s) 


type of escape 


157 


29.4 


axial 


432 


20.2 


axial 


556 


16.1 


axial 


684 


13.7 


radial 


1107 


10.5 


radial 


1108 


11.4 


radial 


1177 


10.7 


radial 


2141 


10.4 


radial 



(35) 



According to table (pq), the atom is so strongly confined axially that in some cases (those with the longer trapping 
times), it leaves the cavity radially, due to the effect of spontaneous emission. In these cases, an additional radial 
cooling mechanism would be helpful as, for instance, that achieved via an additional transverse free-space cooling (see 
pl| ) or via an active feedback on radial motion jljj] . 
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VI. CONCLUSIONS 



We have investigated trapping of single atoms in high finesse optical cavities using both a quasi-resonant cavity 
QED field and an intense classical FORT mode. In particular we have considered the case of doughnut FORT modes, 
and we have compared them to the most common case of a fundamental Gaussian FORT mode (see for example 
the experiment of Ref. [§!). Performing full 3-D numerical simulations of the quasi-classical center-of-mass motion 
of the atom, we have shown that a doughnut FORT mode is more suitable than a fundamental Gaussian FORT 
mode to trap the atom, in the case when the FORT mode shifts both the excited and the ground state down. This 
happens when the FORT mode is red-detuned in such a way that the excited state is relatively closer to resonance 
with a higher-lying excited state than with the ground state. This case, even though more difficult to realize than 
the standard two-level case where the two Stark shifts are opposite, is of particular interest for quantum information 
applications of cavity-QED systems, where it is important to trap the atom at a given position, independently of 
its internal state, so that the quantum manipulation of the internal state can be easily performed. The advantage 
of using a doughnut FORT mode instead of a TEMo.o one is due to the stronger radial confinement, achieved at a 
nonzero distance from the cavity axis. In this way the FORT mode is more suitable to counteract the unavoidable 
radial departure caused by the centrifugal potential and by diffusion, by keeping the atom close to its equilibrium 
position where the cooling provided by the quantized cavity mode is optimal. 
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